New Solutions to the Firing Squad Synchronization Problem 
for Neural and Hyperdag P Systems 



Michael J. Dinneen Yun-Bum Kim 

Radu Nicolescu 

Department of Computer Science, University of Auckland, 
Private Bag 92019, Auckland, New Zealand 
{mjd,radu}@cs . auckland. ac .nz tkim021@aucklanduni . ac ,nz 

We propose two uniform solutions to an open question: the Firing Squad Synchronization Problem 
(FSSP), for hyperdag and symmetric neural P systems, with anonymous cells. Our solutions take 
e c + 5 and 6e c + 7 steps, respectively, where e c is the eccentricity of the commander cell of the dag or 
digraph underlying these P systems. The first and fast solution is based on a novel proposal, which 
dynamically extends P systems with mobile channels. The second solution is substantially longer, 
but is solely based on classical rules and static channels. In contrast to the previous solutions, which 
work for tree-based P systems, our solutions synchronize to any subset of the underlying digraph; 
and do not require membrane polarizations or conditional rules, but require states, as typically used 
in hyperdag and neural P systems. 
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1 Introduction 

The Firing Squad Synchronization Problem (FSSP) (5l|6j|9l[T2l[T4j[T5l is one of the best studied problems 
for cellular automata. The problem involves finding a cellular automaton, such that, after a command is 
given, all the cells, after some finite time, enter a designated firing state simultaneously and for the first 
time. Several variants of FSSP Ifl2l [141 . have been proposed and studied. Studies of these variations 
mainly focus on finding a solution with as few states as possible and possibly running in optimum time. 

There are several applications that require synchronization. We list just three here. At the biological 
level, cell synchronization is a process by which cells at different stages of the cell cycle (division, 
duplication, replication) in a culture are brought to the same phase. There are several biological methods 
used to synchronize cells at specific cell phases [4j|. Once synchronized, monitoring the progression from 
one phase to another allows us to calculate the timing of specific cells' phases. A second example relates 
to operating systems lfl3l . where process synchronization is the coordination of simultaneous threads or 
processes to complete a task without race conditions. Finally, in telecommunication networks Q, we 
often want to synchronize computers to the same time, i.e., primary reference clocks should be used to 
avoid clock offsets. 

The synchronization problem has recently been studied in the framework of P systems. Using tree- 
based P systems, Bernardini et al [2] provided a non-deterministic with time complexity 3h and a de- 
terministic solution with time complexity An + 2h, where h is the height of the tree structure underlying 
the P system and n is the number of membranes of the P system. The deterministic solution requires 
membrane polarization techniques and uses a depth-first-search. 

More recently, Alhazov et al [1] described an improved deterministic algorithm for tree-based P sys- 
tems, that runs in 3h + 3 steps. This solution requires conditional rules (promoters and inhibitors) and 
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combines a breadth-first-search, a broadcast and a convergecast, algorithmic techniques with a high 
potential for parallelism. 

In this paper, we continue the study of FSSP in the framework of P systems, by providing solutions 
for hyperdag P systems [7] and for neural P systems ifTOl with symmetric communication channels. We 
propose deterministic solutions to a variant of FSSP [14], in which there is a single commander, at an 
arbitrary position. We further generalize this problem by synchronizing a subset of cells of the considered 
hyperdag or neural P system. 

These more complex structures pose additional challenges, not considered in the previous FSSP 
papers on tree-based P systems, such as multiple network sources (no single root) and multiple paths 
between cells. Additionally, by allowing an arbitrary position for the commander, we cannot anymore 
take advantage of the sense of direction between adjacent cells; practically, our structures need to be 
treated as undirected graphs. 

Our first solution uses simple rules, but requires dynamical structures. In this paper we propose a 
novel extension, which supports the creation of dynamical structures, by allowing mobile channels. This 
solution works for hP systems and symmetric nP systems; it will also work for tree-based P systems, 
but only if we reconsider them as dag-based P systems, because the resulting structures will be dags, not 
trees. This solution takes e c + 5 steps, where e c is the eccentricity of the commander cell of the under- 
lying dag or digraph. The relative simplicity and the speed of this solution supports our hypothesis that 
basing P systems on dag, instead of tree, structures allows more natural expressions of some fundamental 
distributed algorithms QUI. 

Our second solution is more traditional and does not require dynamical structures, but is substantially 
more complex, combining a breadth-first-search, a broadcast and a convergecast. This solution works 
for tree-based P systems, hP systems and symmetric nP systems and takes 6e c + 7 steps. When restricted 
to P systems, our algorithm takes more steps than Alhazov et al [1J, if the commander is the root node, 
but comparable to this, when the commander is a central node of an unbalanced rooted tree. 

Our two solutions do not require polarizations or conditional rules, but require states, as defined for 
hyperdag and neural P systems 

Section [2] provides background definitions and introduces the families of P systems considered for 
synchronization. Next, in Section |3j we cite the communication models for hyperdag P systems and 
neural P systems, and the transition and rewrite rules available for solving the FSSP. Our two FSSP 
solutions are described in Sections|4]and[5] where we also illustrate the evolution of our FSSP algorithms. 
Finally, we end with some concluding remarks. 

2 Preliminary 

A (binary) relation R over two sets X and Y is a subset of their Cartesian product, RCXxY. For ACX 
andficy, we set R(A) = {y G Y \ 3x€A,(x,y) eR},R^ 1 {B) = {x G X \ 3y£B,(x,y) G R}. 

A digraph (directed graph) G is a pair (X,A), where X is a finite set of elements called nodes (or 
vertices), and A is a binary relation A C X x X, of elements called arcs. A length n — 1 path is a sequence 
of n distinct nodes x\,... ,x n , such that {(jci,X2), . . . , (x„_i,jc„)} Q A. A cycle is a path x\,. . . ,x n , where 
n> 1 and (x n ,xi) GA. A digraph is symmetric if its relation A is symmetric, i.e., (xi,X2) GA<s>(i2,ii) G 
A. By default, all digraphs considered in this paper, and all structures from digraphs (dag, rooted tree, 
see below) will be weakly connected, i.e., each pair of nodes is connected via a chain of arcs, where the 
arc direction is not relevant. 
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A dag (directed acyclic graph) is a digraph (X,A) without cycles. For jcGI, (x) are x's parents, 
A(x) area's children, and A(A~ 1 (x))\{x} are x's siblings. 

A rooted tree is a special case of dag, where each node has exactly one parent, except a distinguished 
node, called root, which has none. 

Throughout this paper, we will use the term graph to denote a symmetric digraph and tree to denote 
a rooted tree. 

For a given tree, dag or digraph, we define e c , the eccentricity of a node c, as the maximum length of 
a shortest path between c and any other reachable node in the corresponding structure. 

For a tree, the set of neighbors of a node x, Neighbor{x), is the union of x's parent and x's children. 
For a dag 8 and node x, we define Neighbor(x) = 8(x) U (x) u5(5"'(x))\{jc}, if we want to include 
the siblings, or, Neighbor(x) = 8 (x) U8~ l (x), otherwise. For a graph G = (X, A), we set Neighbor(x) = 
A(x) = {y | (x,y) G A}. Note that, as defined, Neighbor is always a symmetric relation. 

A special node c of a structure will be designated as the commander. For a given commander c, we 
define the level of a node x, level c (x) G N, as the length of a shortest path between the c and x, over the 
Neighbor relation. 

For a given tree, dag or digraph and commander c, for nodes x and y, if y G Neighbor(x) and 
level c (y) = level c (x) + 1, then x is a predecessor of y and y is successor of jc. Similarly, a node z is 
a /jeer of a node x, if z G Neighbor(x) and level c {z) = level c {x). Note that, for a given node jc, the set of 
peers and the set of successors are disjoint. A node without a successor will be referred to as a terminal. 
We define maxlevel c = max{level c {x) \ x G X} and we note e c = maxlevel c . A level-preserving path from 
c to a node y is a sequence xo, . . . ,X£, such fhatxo = c,x& = y,%i £ Neighbor(xi-i),level c (xi) = i, 1 < i < 
We further define count c {y) as the number of distinct level-preserving paths from c to y. 

The level of a node and number of level-preserving paths to it can be determined by a standard 
breadth-first-search, as shown in Algorithm [T] Intuitively, this algorithm defines a virtual dag based on 
successor relation and, if the original structure is a tree, this algorithm will "reset" the root at another 
node in that tree. 

Algorithm 1 (Determine levels and count level-preserving paths) 

• INPUT: A tree, dag or digraph, with nodes { 1 , . . . , n} and a commander c G {1, . . . ,n}. 

• OUTPUT: The arrays level c \\ and count c [] of shortest distances and number of level-preserving 
paths from c to each node in the structure, over the Neighbor relation. 

array level c [l,. ..,«] = [—1,. .. ,—1]; count c [\,. ..,n] = [0, ... ,0] 
queue Q = () 

Q^c 

level c [c] = 0; count c [c] = 1 
while Q ^ () do 

x^Q 

for each y G Neighborix) do 
if level c [y] = — 1 then 

level c [y\ = level c [x] + 1 
if level c [y] = level c [x] + 1 then 
count c [y] = count c [y] + count c [x] 
return level c 
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Example 1. Figures [TJ [2] and [3] show level c , predecessors, successors, peers and count c , for a tree, a dag 
and a digraph structure, respectively. Small side-arrows indicate the arcs traversed while computing the 
levels, over the induced Neighbor relation, as described in Algorithm [T] 
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Figure 1: Left: a tree (taken from Bernardini et al El), with commander c = 3, ej, =2; Right: table with 
node levels, predecessors, successors, peers and count c 's. 
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Figure 2: Left: a dag with commander c = 6, ^6 = 3 (siblings excluded); Right: table with node levels, 
predecessors, successors, peers and count c 's. 
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Figure 3: Left: a graph with commander c = 1, e\ = 3; Right: table with node levels, predecessors, 
successors, peers and counters. 
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3 P Systems and the Firing Squad Synchronization Problem 

In this section, we briefly recall several fundamental definitions for P systems and describe a P systems 
version of the Firing Squad Synchronization Problem (FSSP). 

For the definitions of tree-based P systems, see Paun [10]. Here we reproduce the basic definitions 
of dag-based hyperdag P systems, from our previous work [7] and digraph-based neural P systems, from 
Paun ifTOl . 

Definition 2 (Hyperdag P systems [7 1) A hyperdag P system (of order n), in short an hP system, is a 
system U h = (0,(Ti, o n ,S,I 0Ut ), where: 

1. O is an ordered finite non-empty alphabet of objects; 

2. Oi,...,a„ are cells, of the form a, = (Qi,Sifi,Wj,o,Pi), 1 < i < n, where: 

• Qi is a finite set (of states), 

• s i,Q ^ Qi i s the initial state, 

• Wifi G O* is the initial multiset of objects, 

• Pj is a finite set of multiset rewrite rules of the form: sx — » s'x'u^ v\w^y so z out , where s, s' G Qi, 
x,x' G O*, Kf G O*, vi G O*, w^ G Ot>, y go G 0* g0 and z out G 0* out , with the restriction that 
Zout = A for all / G { 1 , . . . , n}\I 0Ut . 

3. 8 is a set of dag parent-child arcs on {l,...,n}, i.e., 5 C {l,...,n} x {l,...,n}, representing 
duplex communication channels between cells; 

4- / «f Q {lj •••,"} indicates the output cells, the only cells allowed to send objects to the "environ- 
ment". 



Definition 3 (Neural P systems lHOl ) A neural P system (of order n > 1), in short an nP system, is a 
system n„ = (O, Gi, . . . , a„,5yn, / 0Mf ), where: 

1. is an ordered finite non-empty alphabet of objects; 

2. Oi, ... ,a„ are cells, of the form a,- = (fii^OjWi.Oj^h 1 < i < n, where: 

• Qi is a finite set (of states), 

• Si.o G 2i is the initial state, 

• Wifi G 0* is the initial multiset of objects, 

• is a finite set of multiset rewrite rules of the form: sx — > s'x'ygoZout, where j' G 2;, x,x' G 
0*, y g0 G 0* o and z ou t G 0* HP with the restriction that z 01rt = A for all i G {1, . . . ,n}\{i out }. 

3. syn is a set of digraph arcs on {1, . . . ,n}, i.e., syn C {1, ... ,n} x {1,... ,n}, representing unidirec- 
tional communication channels between cells, known as synapses; 

4- iout £ {lj • • - > n } indicates the output cell, the only cell allowed to send objects to the "environ- 
ment". 

A symmetric nP system, (here) in short, a snP system, is an nP system where the underlying digraph 
syn is symmetric (i.e., a graph). For further definitions describing the evolution of hP and nP systems, 
such as configuration, rewrite modes, transfer modes, transition steps, halting and results, see our previ- 
ous work [ 7]. For all structures, we also utilize the weak policy for applying priorities to rules, as defined 
by Paun O- 
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Remark 4. Most of the P systems considered here (i.e., nP systems, snP systems, hP systems with siblings 
and hP systems without siblings) define a tag go that sends a multiset of objects along the previously 
defined Neighbor relation. Traditional tree -based P systems do not directly provide this facility, however, 
it can be easily provided by the union of out and in\ target indications, that represent sending "to parent" 
and "to all children", respectively. That is, (w,go) = (w,out)(w,in\). 

Definition 5 (FSSP for P systems with states — informal definition) We are given a P, hP, snP or 

nP system with n cells, {Oi , . . . , a,,}, where all cells have the same states set and same rules set. Two 
states are distinguished: an initial state so and a. firing state s$. We select an arbitrary commander cell 
(7 C and an arbitrary subset of squad cells, F C {oj,. . . ,(T B } (possibly the whole set), that we wish to 
synchronize; the commander itself may or may not be part of the firing squad. At startup, all cells start 
in the initial state s$; the commander and the squad cells may contain specific objects, but all other cells 
are empty. Initially, all cells, except the commander, are idle, and will remain idle until they receive a 
message. The commander sends one or more orders, to one or more of its neighbors, to start and control 
the synchronization process. Idle cells may become active upon receiving a first message. Notifications 
may be relayed to all cells, as necessary. Eventually, all cells in the squad set F will enter the designated 
firing state s$ , simultaneously and for the first time. At that time, all the other cells have reached a 
different state, typically sq or s\, without ever passing through the firing state s$. Optionally, at that time, 
all cells should be empty. 

In this paper, we propose two new deterministic FSSP solutions, that are described in the next two 
sections. Our two solutions do not require polarities or conditional rules, but require priorities and states. 
Both hP systems and snP systems already have states, by definition. However, it seems that traditional 
tree-based P systems have not used states so far, or not much. 

4 FSSP — Dynamic Structures via Mobile Channels 

In this section, we further refine our solution given in an earlier paper HQ. A natural solution is possible 
when we are allowed to extend the cell structure of the given hP or snP system. We achieve this by 
supporting mobile channels. The endpoints of our mobile channels appear in the rules like all other 
objects and are subject to usual rewriting and transfer rules. The end result is a channel that grows 
step-by-step, not unlike a nerve which extends in a growing or regrowing tissue. 

We first extend the original hP or nP system by an external cell, which will be called the sergeant. 
Next, this sergeant will send a self -replicating mobile endpoint, that will be repeatedly broadcasted, until 
all cells are reached. A mobile endpoint will leave a fixed endpoint in a squad cell and disappear without 
trace from the other cells. In the end, the structure will be extended with new channels which will link 
the sergeant with all squad cells. Finally, when there are no more structural changes, the sergeant, will 
send a firing command to all squad cells, prompting these cells to enter the firing state, all at the same 
time. 

Our algorithm uses the following special objects, which can in principle be rewritten and transferred 
as all other traditional objects, but, at the same time, are also endpoints for dynamically created channels: 

• a is here the fixed endpoint of all dynamically created channels (here we use only one a) 

• 6 is a mobile endpoint of a dynamically created channel (here this symbol is further processed by 
rewriting and transfer rules) 

• CO is a new fixed endpoint of a dynamically created channel (here this symbol will remain fixed) 
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Briefly, the initial structure is dynamically extended by all arcs (a a ,a t ), where o a is the cell that contains 
a and a t is a cell that contains 6 or co. 

The following algorithm assumes that the first step has already been completed, i.e., the sergeant was 
already created by one of the existing cell creation or division rules, already available for P systems. 



Algorithm 2 (FSSP — Dynamic structures via mobile channels) 



Precondition: An hP or nP system, with n cells G\, . . . , a n , a commander cell o c and a set of squad cells 
F to be synchronized. Additionally, we already have a sergeant cell, G n+ \, linked to the commander by 
one arc, (a n+ i,(7 c ), for hP systems, or by two arcs, (o n+ \,o c ), (<t c ,<t„ + i), for snP systems. 

All cells start in the state so and have the same rules. The state is the firing state. Initially, the 
sergeant o„+\ is marked by one object a, and each squad cell is marked by one object / (this can include 
the commander o c or the sergeant o n+ \, or both); all other cells have no objects. 

Postcondition: All cells in the set F enter state sq, simultaneously and for the first time, after e c + 5 
steps, where e c is the commander's eccentricity in the underlying graph. All other cells enter state si, 
without ever passing through state ^ . 



Rules (rules are applied under the weak interpretation of priorities, in the rewrite mode a 
transfer mode j8 = re pi): 
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Example 6. Figure [4] and Table [T] illustrate this algorithm for an hP system based on the dag of Figure [T] 
Here, the commander cell is 03, the squad set is F = {o\, . . . , 05} and this system's structure has already 
been extended by the sergeant cell as and the arc ((Tg, 03). The mobile channels are represented by dotted 
arrows. 
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Figure 4: Running Algorithm [2] on the hP system of Example[6] 
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Table 1: Traces for the hP system of Example |6] 
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5 FSSP— Static Structures and Rules 

Here we consider a second scenario, where we are allowed to modify the rules of the given hP or nP sys- 
tem, but not its original structure. A brief description of this solution follows. The commander intends to 
send an order to all cells in the set F, which will prompt them to synchronize by entering the designated 
firing state. However, in general, the commander does not have direct communication channels with all 
the cells. In this case, the process of sending a command to the destination cell will cause delays (some 
steps), as the command is relayed through intermediate cells. Hence, to ensure all firing squad cells enter 
the firing state simultaneously, each firing squad cell determines the number of steps it needs to wait 
before entering the firing state. 

As in our earlier paper [8], cells have no built-in knowledge of the network topology. Additionally, 
cells are anonymous, i.e., not identified by cell IDs, and not implicitly named by membrane polarization 
techniques. The cells are initially empty, except the commander, which is initially marked by one a, and 
the squad cells, which are initially marked by one / each. All cells start with the same set of rules, which 
are applied in the max rewrite mode, using weak priorities, and the repl transfer mode. In the proofs, 
all rules that are concurrently applied will be grouped together within parentheses; e.g., (x,y),z indicates 
two steps, first rules x and y, concurrently executed, followed by rule z. 

Each cell independently progresses fhrough/owr phases, called FSSP-I, FSSP- II, FSSP-III and FSSP- 
IV, which are detailed in Algorithms [3] [4] [5] and [6j respectively. An overview of these four phases is as 
follows: 

• Phase FSSP-I is a broadcast from the commander, that follows the virtual dag defined by level c . 
This phase starts in state so and ends in state S2- Also, the commander starts a counter, which, at 
the end of Phase FSSP-II, will determine its eccentricity. 

• Phase FSSP-II is a subsequent convergecast from terminal cells, that follows the same virtual dag. 
This phase starts in state S2 and ends when the commander enters state s^,. At the end of this phase, 
the commander's counter determines its eccentricity. 

• Phase FSSP-III is a second broadcast, initiated from the commander, that follows the same virtual 
dag. This phase starts in state S(, and ends in state s&. The commander sends out its eccentricity, 
which is successively decremented at each level. 

• Phase FSSP-IV is a timing (countdown) for entering the firing state. This phase starts in state s& 
and continues with a countdown, until squad cells simultaneously enter the firing state sg, and all 
other cells enter state sq. 
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The statechart in Figure [5] illustrates the combined flow of these four phases. The nodes represent 
the states of the hP or nP system and the arcs are labelled with numbers of the rules that match the cor- 
responding transitions. The rest of this section describes these four phases, proving their correctness and 



time complexities. A sample run of our algorithm will follow at the end of this section, in Example 12 



2.1,2.7 



" 6) 6.1,6.2"® 1 7.2,7.3,' ,| S^ 8.2,8.4 




Figure 5: Statechart view of the combined FSSP algorithm phases. 



FSSP: The initial configuration 



r = {oi, . . . , o n }, n> 1, is the set of all cells, o c is the commander, and the firing squad is F C T; 



O = {a,b,c,d,e,f,g,h,k,l,p,q}; 



• Qi = {*0) s\, *2) ^3) ^5) ^7) *8) J 9}» for i £ {1, . . • ,«}, which is "allocated" to four phases as fol- 
lows: FSSP-I contains rules for states {so^i}; FSSP-II contains rules for states {52,53,54,55,56}; 
FSSP-III contains rules for states {56,57}; FSSP-IV contains rules for states {58,59}; 



50 = 59 is the firing state; 



slo = so, fori e{l,...,n}; 



• w c .Q = {a}, if o c F, or {a,f}, otherwise; w,,o = {/} for all a,- G F\ o c ; Wifi = 0, for all a,- G 
r\(FU{a f }); 



• The following rules are applied under the weak interpretation of priorities, in the rewrite mode 
a = min and transfer mode j8 =repl: 
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4. For state 54: 

1) s^cd — > 54 
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5. For state 55: 

1) 5 5 &^5 5 

2) 55a — > s^ac g0 

3) 5 5 /jp -> 5 5 p 

4) s 5 pq^s 5 

5) s 5 p^s 5 kp 

6) 55/ -> 5 5 //z g0 

7) * 5 / -> s 6 q g0 



Algorithm 3 (FSSP-I: First broadcast from the commander) 



6. For state ^ 

1) S(,ae — ► 57«fc 

2) 5 6 e -> 57 

3) 5 6 c -> 5 6 

4) 5 6 g^5 6 

5) 5 6 /j -> 5 6 

6) S 6 /7 -> 5 6 

7) 5 6 g^5 6 

7. For state 57: 

1) Sqk — > 57 

2) 57a — * s%a 

3) 5 7 e -> 5 8 

8. For state 5§: 

1) 58^ — > 58a 

2) 58a/ — > 59 

3) 5 8 a ->• 5 

4) 58a — > 59 



Precondition: The initial configuration as specified earlier. 
Postcondition: 

• The end state is 52. 

• A cell a, has 

o count c (i) copies of a and count c (i) copies of k; 

o u copies of /, where u is the total number of a's in a,'s peers; 

o v copies of d, where v is the total number of a's in cr,'s successors; 

o two copies of e, if a, = G c ; 

o one copy of /, if a, G F. 



Proof. This phase of the algorithm is a broadcast that follows the virtual dag created by the levels 
determined by Algorithm [T] 

Consider a cell a,. By induction: 

• At step level c (i), a, (except the commander) receives a total of count c (i) copies of d from its 
predecessors. 
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• At step level c (i) + 1, a, broadcasts count c (i) copies of d to each of its neighbors and transits to 
state s\. At the same time, a, accumulates one local copy of a for each sent d, for a total count of 
count c (i) of a's. Also, a, receives u copies of d, similarly sent by its peers, where u is equal to the 
total number of a's similarly accumulated, at the same time step, by a,'s peers. 

• At step level c (i) + 2, a, receives v copies of d, sent back by its successors; and transits to state S2, 
where v is equal to the total number of a's created, at the same time step, by O/'s successors; 

The commander, by initially having one a, creates two copies of e. Finally, the rules associated with 
this phase do not change the number of /'s, thus, each cell in the firing squad still ends with one /. □ 

Corollary 7 (FSSP-I: Number of steps). For each cell Ou the phase FSSP-I takes level c (i) + 2 steps. 
Proof. As indicated in the proof of the Algorithm|3} the total number of steps is level c (i) + 2. □ 



Algorithm 4 (FSSP-II: Convergecasts from terminal nodes) 



Precondition: As described in the postcondition of Algorithm [3] 
Postcondition: 

• This phase ends when the commander enters state s^,. 

• A cell a, has 

o count c {i) copies of a; 

o e c + 2 copies of e, if a, = a c ; 

o one copy of /, if a, G F. 

Proof. Briefly, this phase of the algorithm is a convergecast of c's, starting from terminal cells, and 
further relayed up, on the virtual dag, until the commander is reached. 

For the purpose of this phase, the non-commander cells can be organized in the following three 
groups: TC cells = terminal cells; NTC-NTP cells = non-terminal cells without non-terminal peers (i.e., 
cells without peers or cells with terminal peers only); NTC+NTP cells = non-terminal cells with non- 
terminal peers (these cells may also have terminal peers). 

During this phase, these cells will make transitions between the following three conceptual stages: 
WCS = waiting for convergecasts from successors (state S4); RTC = ready to convergecast (state S5); HC 
= have convergecasted (state s 6 ). Specifically, the following transitions will be made: the TC cells will 
transit immediately from the WCS stage to the HC stage; the NTC-NTP cells will linger in the WCS 
stage until they receive convergecasts from all their successors, after which they will transit directly to 
the HC stage; the NTC+NTP cells will linger in the WCS stage until they receive convergecasts from all 
their successors, subsequently they will linger in the RTC stage until all their non-terminal peers reach 
the RTC stage as well, after which they will transit to the HC stage. 

During this process, cells will exchange c-notifications, which are messages consisting of number of 
c's and /z-notifications, which are messages consisting of number of h's. The actual numbers depend on 
network topology and take into account the multiple paths that appear in the virtual dag. 

The c-notincation broadcasted by cell a, consists of count c (i) copies of c and is only sent once when 
a, transits into the HC stage. 

The /j-notification broadcasted by cell a, consists of u copies of h, where u is the number defined 
in the precondition. This notification is sent repeatedly, while a, remains in the RTC stage, until a, 
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transits into the HC stage. The /j-notirications synchronize the non-terminal peers that cannot transit to 
the HC stage until all of them have reached the RTC stage. This avoids the potential confusion that could 
otherwise arise when a non-terminal cell receives an "ambiguous" c-notification, i.e., a c-notification that 
could come both from a successor or from a non-terminal peer. 

Without loss of generality, we illustrate our solution on the dag from Figure [6j This figure shows a 
typical sub-dag of the virtual dag created by Algorithm [T] where cells G\, 02, 03, C4 and 05 are at the 
same level, and the horizontal lines indicate peer relations. 




OOOO 

Figure 6: A typical sub-dag of the virtual dag. 

Cell (7i is a TC cell and will transit from stage WCS to stage HC, immediately after detecting that it 
has no successors (there is no need for any synchronization with its peer G2). 

Cell G5 is a NTC-NTP cell and will transit from stage WCS to stage HC, after receiving c-notifications 
from all its successors. 

Cells 02, 03 and 04 are NTC+NTP cells. Each of these cells will linger in stage WCS until it receives 
c-notifications from all its successors, when it will enter stage RTC. Cells 02 and 03 are peers, therefore 
none of them will be allowed to transit to stage HC, until both of them have reached the RTC stage. 
Similarly, cells 03 and 04 are peers, therefore none of them will be allowed to transit to stage HC, until 
both of them have reached the RTC stage. Assume that cells 02, 03 and 04 will reach stage RTC in this 
order. Then, cell 02 will wait in stage RTC until 03 also reaches the same stage. When this eventuates, 
(T2 will transit to stage HC, while 03 will still linger in stage RTC until 04 reaches the same stage. When 
this eventuates, both 03 and 04 will transit at the same time to stage HC. 

A TC cell (7,- enters this phase level c (i) steps after the commander, idles one step in state S2, then 
starts its role in the convergecast, by broadcasting count c (i) copies of c to its predecessors and peers (it 
does not have successors) and transits to state This cell further idles in state S6 until it receives e's 
from its predecessors. The convergecast takes four steps at each level. 

The total run-time is dominated by e c , the length of the longest level-preserving path from comman- 
der. Therefore, the convergecast wave will complete at commander after e c +4e c — 2 = 5e c — 2 steps after 
the commander starts this phase. When the commander receives the convergecast from all its successors, 
it takes two steps to transit to state s$. Therefore, the commander enters state s&, 5e c steps after it starts 
this phase. □ 

Corollary 8 (FSSP-II: Number of steps). For each cell a„ the phase FSSP-II takes 5e c — level c (i) steps. 



Proof. As indicated in the proof of Algorithm|4j this phase takes 5e c steps. 



□ 
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Algorithm 5 (FSSP-III: Second broadcast from the commander) 

Precondition: As described in the postcondition of Algorithm [4] 
Postcondition: 

• The end state is s% . 

• A cell a, has 

o count c {i) copies of a; 

o {e c + 1 — level c {i))count c {i) copies of b; 

o one copy of /, if a* e F. 



Proof. In this phase, commander starts its second broadcast, by sending e c + \ copies of e's to all its 
successors. By induction on level, a cell a, receives a total of [e c + 2 — level c {i))count c {i) copies of e's 
from its predecessors, reduces this count by count c (i) (i.e., the count of a's), forwards the remaining (e c + 
1 — lev el c (i)) count c (i) copies of e's to all its successors and creates for itself [e c + 1 — lev el c {i)) count c {i) 
copies of b's. A more detailed description will be given in the final version. 

All rules of this phase do not change the number of a's or the number of f's; therefore, the corre- 
sponding postcondition holds. rj 

Corollary 9 (FSSP-III: Number of steps). For each cell Ou the phase FSSP-III takes level c (i) + 3 steps. 

Proof. As indicated in the proof of Algorithm|5} this phase takes level c (i) + 3 steps. □ 

Algorithm 6 (FSSP-IV: Timing for entering the firing state) 

Precondition: As described in the postcondition of Algorithm [5] 
Postcondition: 

• The end state is sg for cells in the firing squad, or so, otherwise. 

• Each cell is empty. 



Proof. As long as b's are present, a cell a,- performs a transition step that decreases the number of b's by 
count c {i) (i.e., the number of a's). This step will be repeated [e c + 1 — level c {i)) times, as given by the 
initial ratio between the number of b's, [e c + 1 — lev el c {i)) count c {i), and the number of a's, count c (i). 
This is the delay every cell needs to wait, before entering either the firing state sg or the initial state so. 

Finally, in the last step, cell (7/ enters sg, if a, has one f, or so> otherwise. At the same time, all 
existing objects are removed. □ 

Corollary 10 (FSSP-IV: Number of steps). For each cell a,-, the phase FSSP-IV takes e c + 2 — level c {i) 
steps. 



Proof. As indicated in the proof of Algorithm |6j this phase takes (e c + 1 — level c {i)) + 1 = e c + 2 — 
level c (i). □ 
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Theorem 11. For each cell a,, the combined running time of the four phases Algorithm^ [4] [5]a7iG?[6]zs 

6e c +7, where e c is the eccentricity of the commander o c . 

Proof. The result is obtained by summing the individual running times of the four phases, as given by 
Corollaries |7J [8} [9] and [To} (level c (i) + 2) + (5e c - level c (i)) + (level c (i) + 3) + (e c + 2 - level c (i)) = 
6e c + l. □ 

Example 12. We present traces of the FSSP algorithm for the hP system given in Figure [2] in Table [2j 
where the cells are ordered according to their levels and the starting states of phases FSSP-II, FSSP-III 
and FSSP-IV are highlighted. 



Table 2: The FSSP trace on the dag of Figure[2| where c = 6, <?6 = 3, F = {o\ , 04, 05, 07, 09, Cio}. 
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6 Conclusion 

We have presented two new algorithms for the Firing Squad Synchronization Problem that operate on 
several families of P systems. Out of the box, both algorithms work for hyperdag P systems and sym- 
metric neural P systems. The first algorithm is based on dynamic structures and highlights the merits of 
dags as underlying structures for P systems. To support the required dynamic structures, we propose an 
extended interpretation of P systems which allows mobile channels, a solution which we believe is fully 
compatible with the existing P systems rules. 

The second algorithm, which is more complex, is applicable to P systems with static membrane 
topologies and is uniformly defined in terms of a structural Neighbor relation. These two algorithms do 
not require naming facilities, such as cell IDs or cell polarization and handle a generalized version of 
the FSSP, where the commander can assume an arbitrary position and only a specified subset of the cells 
needs to be synchronized. 

The work started in this paper leaves open several interesting problems. Can we find simpler and 
more efficient solutions for hP systems based on single-sourced dags? Can we find simpler and more 
efficient solutions for hP or snP systems using named cells (unique cell IDs)? Can we find a solution 
for arbitrary strongly-connected (non necessarily symmetric) nP systems? What is relation between the 
mobile channels, which we have here proposed for P systems, and the support for mobile channels in the 
71-calculus? 
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